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We investigate the time-of-night variation of solar neutrino rate which will be of relevance to Super- 
Kamioka and Sudbury neutrino detectors in the framework of oscillations among the three flavors. 
An analytical method of computing the regeneration in the earth is presented. If day-night effect 
is seen, we show how the study of the time-of-night variation will allow the determination of the 
neutrino parameters. 

PACS number: 14.60.Pq, 95.30.Cq, 96.40. Tv, 96.60.Jw 



It has been known for quite some time that an asym- 
metry between night rate and day rate for the real time 
solar neutrino detectors is an unambiguous signal for neu- 
trino mixing and oscillations. Conversely the absence of 
such an effect can put constraints on the neutrino masses 
and mixing angles. Although no day-night asymmetry 
outside the error-bars was seen at the Kamioka detector 
|y,|) the high statistics detectors like Super-Kamioka j| , 
SNO ||]and Borexino ||] will be much more effective in 
investigating this effect. If there is a day-night asymme- 
try, then the profile of the asymmetry as a function of 
the time during night is a very sensitive function of the 
neutrino parameters. The counting rates in these detec- 
tors are expected to be high enough for the study of this 
time of night variation. The aim of this note is to fo- 
cus attention on this aspect of the day-night effect Q , in 
view of the fact that Super-Kamioka has already started 
functioning and SNO is expected to do so soon. 

The neutrino samples different amounts of matter in 
the earth during a single night and also during a year. 
The distance d travelled by the neutrino inside the earth 
during night, as a function of time t, is given by 

2nt 
d = 2R(smcf)i sin S + cos 4>i cos<5cos(— — )) , (1) 

Td 



2nt 
sin 8 = sin 23.5° sin( ), 
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R is the radius of the earth, (f>i is the latitude of the loca- 
tion of the detector, Tjj is the length of the day, Ty is the 
length of the year and zero of t is chosen at midnight on 
autumnal equinox. Thus assuming that the neutrino pa- 
rameters are in a suitable range, neutrino data collected 
as a function of t contain an enormous amount of infor- 
mation on neutrino oscillations, which in principle can 
be analyzed to yield the neutrino parameters. 

The time variation of x = ( — ) during the night and 
n ^-R 
year are illustrated in Fig.pla. If the data collected during 

successive nights are accumulated at the corresponding 

x— bins, the calculation of neutrino rates per unit bin will 

require the function f(x) defined as the time duration 

per unit interval of a; H. f(x) for different locations are 

plotted in Figjjlb, which shows the relative merits of the 
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FIG. 1. (a) The fractional distance x travelled by the neu- 
trino inside the earth during night is plotted against time 
for the latitude (36.4°) of the Super-Kamioka detector. The 
lower figure gives the envelope of the 365 maxima during the 
year . As examples of the actual curves, those for a few nights 
during three specific seasons of the year are shown in the up- 
per figure. (b) The function f(x) in hours per unit x is 
plotted for various latitudes: Super-Kamioka (36.4°), Borex- 
ino (42.45°), SNO (46.5°) , equator (0.0°) , pole (90°). 



detectors for exposure to regions of x. 

We now describe an analytical method of calculating 
the neutrino regeneration effect in the earth. Let a neu- 
trino of flavor a be produced at time t = to in the core 
of the sun. Its state vector is 



i^mmhkhE^ 



c ^?) 



(3) 



where \vf) are the mass eigenstates with mass eigenval- 
ues nf and U^i are the elements of the mixing matrix 
in the core of the sun. We use Greek index a to denote 
the flavors e, fj,,r and Latin index i to denote the mass 
eigenstates i = 1,2,3. The neutrino propagates in the sun 
adiabatically upto tn (the resonance point), makes nona- 
diabatic transitions at tji, propagates adiabatically upto 
t\ (the edge of the sun) and propagates as a free parti- 
cle upto £2 when it makes a nonadiabatic transition as it 
enters the earth. If it propagates adiabatically upto £3 
(we shall soon correct for nonadiabatic jumps during this 
propagation), the probability amplitude for detecting a 
neutrino of flavor j3 at £3 is 



k,j,i 
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s E {t)dt 



+ e 3 (t 2 -h)+ J e s 3 {t)dt + \ ef(t)dt ) } (4) 



£j(= E + p 2 /2E) ,ef(t) and ef(t) are the energy eigen- 
values in vacuum, sun and earth respectively, M^ and 
M E t are the probability amplitudes for nonadiabatic tran- 
sition i — > j inside the sun and at the vacuum -earth 
boundary respectively. The latter is due to the abrupt 
change in density when the neutrino enters the earth and 
is given by 



M% = {^\v 3 ) = 5>f k)Kh-> = J2 U ?kU* aj (5) 



where U and U E are the mixing matrices in vac- 
uum and earth respectively. Averaging the probability 
K^pl^aih))] 2 over tR results in the desired incoherent 
mixture of mass eigenstates of neutrinos reaching the 
surface of the earth. Calling this averaged probability 
as P^p ( the probabilty for a neutrino produced in the 
sun as v a to be detected as vp in the earth at night), we 
can write the result as 



i 
^• = £|M|| 2 |t/£| 2 , (7) 

i 

p % = E UftUfcMOJiMfiexp (-2t* kfc ,) , (8) 

(9) 



fc.fc' 



**v = 2j (4(t)-4>(t))dt 



For the daytime, put i 3 = t 2 so that F% becomes \Upj\ 2 
and so eq.(y) reduces to the usual |10| , [llfl transition prob- 
ability in the day: 
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It is important to note that the factorization of prob- 
abilities seen in eq.(|l) is valid only for mass-eigenstates 
in the intermediate state. An equivalent statement of 
this result is that the density matrix is diagonal only in 
the mass-eigenstate representation and not in the flavor 
representation. 

We next show how to take into account nonadiabatic 
jumps during the propagation inside the earth. Consider 
v propagation through a series of slabs of matter, density 
varying inside each slab smoothly but changing abruptly 
at the junction between adjacent slabs. The state vector 
of the neutrino at the end of the n th slab \n) is related 
to that at the end of the (n — l) th slab \n — 1) by \n) — 
F^ n 'M^ n '\n — 1) where M^ describes the nonadiabatic 



jump occuring at the junction between the (n — 1) and 
n th slabs while F^ n > describes the adiabatic propagation 
in the n th slab. They are given by 

M% ] = {vt i] \vf- l) ) = (t/^V"- 1 ^. , (11) 
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Sijexp I -i / Ei(t)dt J , 



(12) 



where the indices (n) and (n — 1) occuring on v and U 
refer respectively to the n th and (n — l) th slabs at the 
junction between these slabs. Also note that M^ is the 
same as M E defined in eq(B). Defining the density matrix 
at the end of the n th slab as p( n ' — \n){n\ , we have the 
recursion formula 



p(") _ p{n) j^{n) ln-X) j^-{n)\ p{n)\ ^ 



(13) 



Starting with p(°' = \Vj){i/j\ (i.c Vj entering the earth), 
we can calculate p( N 1 at the end of the N th slab using 
(fut) . The probability of observing vp at the end of the 
N^ 1 slab is 



i/3 



= {vp\pW\vp) = {UWpWU^)pp. (14) 



This formula (which reduces to eq(@) for N — 1) can be 
used for the earth modelled as consisting of (iV+l)/2 con- 
centric shells, with the density varying gradually within 
each shell. We shall present numerical results for N = 3 
(mantle and core) |lSq| . However for x < 0.84, neutri- 
nos pass only through the mantle and so N = 1. Accu- 
racy achieved with this model is adequate for the present 
purposes, but the formalism allows one to improve the 
accuracy to any desired level, by adding more shells. 

Apart from the nonadiabatic jumps occuring at the 
density-discontinuities, such jumps can occur also at any 
MSW resonance in the earth. The formalism presented 
above is capable of handling this. One simply replaces 
eq.(|Tl]) for MW for that transition by an appropriate 
Landau-Zener formula ||l3| . 

We parametrize the mixing matrix U in vacuum as 
U = U 23 (ip)U 13 ((f>)U 12 (Lu) where U lj {8 l0 ) is the two fla- 
vor mixing matrix between the ith and the jth mass eigen- 
states with the mixing angle 0y, neglecting CP violation. 
In the solar neutrino problem ip drops out |0,H. The 
mass differences in vacuum are defined as S 2 i = p 2 — p\ 
and J31 = pi — p\. It has been shown [ pTpq | that the 
simultaneous solution of both the solar and the atmo- 
spheric neutrino problems requires the mass hierarchy 
831 3> #21 and under this condition J31 also drops out. 
The rediagonalization of the mass matrix in the presence 
of matter (in the sun or earth) under the hierarchy con- 
dition leads to the following results [[Llfl 



tan 2lj„-i = 



621 sin 2uj 



J21 cos 2oj — A cos 2 <f> ' 
sin 4> m — sin , 
^21 = ^21 cos 2(w — u) m ) — A cos 2 1 



1 cos 2u) r , 



(15) 

(16) 
(17) 
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FIG. 2. A as a function of x, the fractional distance travelled by the neutrino, for stated values of to, S21 (in eV 2 ) and 



where A is the Wolfenstein term A = 2\f2 G F N e E (N e 
is the number density of electrons and E is the neutrino 
energy) . We note that 831 ^> A, for A evaluated at any 
point in the sun or the earth. In eqs ( |l5| ) - ("J]) , the "to" 
stands for matter and in using these equations, one must 
use the appropriate density of matter that is required at 
the various points along the trajectory of the neutrino. 
All the probabilities P^g, P^ and PR satisfy the nor- 



malization conditions, as for instance, J2jPja = 1- for 



three flavors, use of these conditions allows us to express 
P e N e in terms of Pg, Pf e , P 2 E e and P e % as 



cos 2 (f)(sm 2 usP le 



-^((3 cos 2 cos 2 w- l)Pf e - (3cos 2 
— cos 2 4> cos 2 to)] I (cos 2 4> cos 2w) . 



cos 2 uPE) 



)sin ui 



l)Pl E e 

(18) 



This is a general formula that goes over to the one given 
in Ref. |lq| for two flavors. Simplifications arise from 
the mass heirarchy condition. First, My is nonzero for 
i,j = 1,2 only and hence we can replace P^ 3 in eqjlq) 
by sin 2 cf>c- | Af^| 2 is taken to be the modified Landau- 
Zener jump probability for an exponentially varying solar 
density p3[ . Further, M^ also are reduced to 2 x 
2 matrices. In fact eq(pd|) gives M^ = U 12 (9) where 
U 12 (9) is the 2-flavor mixing matrix, with mixing angle 



6 = u) n — uj„-i- As a result, we get the simple formulae 
from cq(||), valid for x < 0.84 : 

Pf e = cos 2 0[cos 2 uj e - sin 2u>e sin 2(ui E - ui) sin 2 $i 2 ] (19) 
P^ e = cos 2 0[sin 2 u E + sin 2ui E sin 2(u E - u>) sin 2 $i 2 ] (20) 

where uj e is the mixing angle, just below the surface of 
the earth. 

The neutrino detection rates for a Super-Kamioka type 
of detector is given by 



R 



baP ee dE + ^ J -oil. /',, xl E (21) 



where the second term is the neutral current contribu- 
tion and (j){E) is the solar neutrino flux as a function 
of the neutrino energy E and cr(E) is the cross section 
from neutrino electron scattering and we integrate from 
5MeV onwards. The cross section is taken from |17| ] 
and the flux from |lq |. The rates for the night and 
day .Rat and Rd are calculated using P^ e and P® re- 
spectively. We define the day-night asymmetry ratio as 
A = (Rjy — Rd)/{Rn + Rd)- We can multiply Rn and 
Rd by the function f{x) displayed in Fig.|l| to get the 
rates per unit interval in x. Note that f{x) cancels in 
the asymmetry ratio calculated theoretically. However, 



the experimentalists have to weight the day rate Rd with 
f(x) before comparing their data with our theoretical 
curves. 

In Fig.g and g we have plotted A as a function of x, the 
fractional distance travelled by the neutrino inside the 
earth for various values of the neutrino parameters 621,00 
and <p. Different values of these parameters have distin- 
guishable characteristics. Some gross features which may 
enable us to specify their approximate domains are the 
following: 

• For small angle u there is a gradual increase of the 
asymmetry with x, whereas for large u> the oscil- 
lations in x start showing up. For x < 0.84 (i.e. 
trajectories through mantle only) there is a very 
clear discrimination between the small w and large 
uj, irrespective of 621 and <f>. 

• As increases, the asymmetry at any x decreases. 

• The amplitude of the oscillatory pattern is largest 
for small £21 and decreases steadily as £21 increases. 

• For small u) and large 82 1, asymmetry is appreciable 
only in the core and is a sensitive function of <52i ■ 

For non-zero <fi a fraction of the solar neutrinos come 
out of the sun as v T and these v T cannot reconvert back 
to v e inside earth because <f> is not affected by matter 
(see eqn.(16)). So a non zero <fr dilutes the asymmetry. 
Our numerical results include the effect of any adiabatic 
MSW resonances that may occur inside the earth. For 
<f> = 0, as pointed out recently [|]|J, MSW resonances 
do occur in the earth's core. However, for large <$> they 
disappear and this is another reason for the regeneration 
efect to be smaller for large <fi, in the core. 

In Fig. I) we have chosen a few parameter sets for 
which A is very small (< 0.15) since they are possible so- 
lutions to the solar and atmospheric neutrino problems 
pl] , p~5|JT9| . But we cannot rigorously exclude other val- 
ues of the neutrino parameters at the present stage of 
knowledge. Day-night effect must be studied in an unbi- 
ased manner, especially because the ratio A is relatively 
independent of the uncertainties of the solar models. 

After this work was completed, we came to know 
of a related work: E. Lisi and D. Montanino, hcp- 
ph/ j)702343| . They have also stressed the importance of 
an analytic approach, however their method is different 
from ours and they have ignored the third flavor. 

We thank M. C. Sinha, M. V. N. Murthy and S. Uma 
Sankar for discussions. 
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FIG. 3. A as a function of x, for (a) 621 = 3 x 10" 5 eV 2 
w = 28°, = 0; (b) 821 = 6 x 1(T 6 , uj = 2.5°, = 0; (c) 8 2 
= 2.5 x 10~ 5 , w = 20°, = 32°; (d) 621 = 6.3 x 10~ 5 , uj = 
18°, = 45°; (e) 8 21 = 2.5 x 10 -6 , u = 6°, = 50°. 



[2] 
[3] 

[4] 

[5] 



[6] 
[7] 

[8] 

[9] 
[10] 

[11] 



[12] 
[13] 
[14] 
[15] 
[16] 

[17] 
[18] 
[19] 



Kamiokande Collaboration: Y.Fukuda et aZ,Phys. Rev. 
Lett. 77, 1683 (1996). 

Y.Totsuka, Nucl. Phys.(Proc. Suppl.) B48, 547 (1996); 
A. Suzuki, in Physics and Astrophysics of Neutrinos, 
edited by M. Fukugita and A. Suzuki ( Springer- Verlag, 
Tokyo, 1994), page 414. 

G.T.Ewan et a/,SNO Collabaration, "Sudbury Neutrino 
Observatory Proposal", Report No. SNO-87-12, 1987 
(unpublished) . 

R.S. Raghavan, in Proceedings of the 25th International 
Conference on High Energy Physics, Singapore, 1990, 
edited by K.K. Phua and Y. Yamaguchi (World Scien- 
tific, Singapore, (1991), p. 698. 
For earlier work, see M-ity and references therein. 
M. L. Cherry and K .Lande, Phys. Rev, D36, 3571 (1987) 
A. J. Baltz and J. Weneser, Phys. Rev, D35, 528 (1987); 
i&«Z,D50,5971 (1994); ibid, D51, 3960 (1994) 



J. M. Gelb, W. Kwong, and S. P. Rosen, tiep-ph/9612332 



Stephen Parke, Phys. Rev. Lett. 57, 1275 (1986). 

M. Narayan, M. V. N. Murthy, G. Rajasekaran and S. 

Uma Sankar, Phys. Rev. D53, 2809 (1996), ibid. D56, 

437 (1997). 

We have taken the densities from F. D. Stacey, Physics 

Of the Earth, John Wiley & Sons, Inc., New York. (1977). 

T. K. Kuo and J. Pantaleone, Rev. Mod. Phys. 61, 937 

(1989). 

A. S. Joshipura and M. V. N. Murthy, Phys. Rev. D37, 

1374 (1988) 

G.L.Fogli, E. Lisi and D. Montanino, Phys. Rev. D49, 

3626 (1994);i6?rf. D54, 2048 (1996). 

S.P.Mikheyev and AYu. Smirnov, in New and Exotic 

Phenomena, Proceedings of the Moriond Workshop, Les 

Arc Savoie, France, 1987, edited by O.Fackler and J. 

Tran Thanh Van (Editions Frontie res, Gif-Sur-Yvette, 

France, 1987), p. 405. 

J.Bahcall, Neutrino Astrophysics, Cambridge University 

Press,New York, (1989) 

J. N. Bahcall and M. H. Pinsonneault, Rev. Mod.Phys. 

64,885 (1992). 

N. Hata and P. Langacker, Phys. Rev. D52, 420 (1995). 



[1] Kamiokande Collaboration: K.S.Hirata et al,Phys. Rev. 
D 44, 2241 (1991). 



